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1 Prologue 



The AdS/CFT correspondence relates string states and local gauge-invariant operators 
of a dual quantum field theory 0. The energies of the string states correspond to the 
eigenvalues (the so-called anomalous dimensions) of the mixing matrix of gauge field 
theory operators. However, even for the well understood case of the M = 4 super 
Yang-Mills (SYM) theory, testing of such a correspondence has revealed to be a rather 
difficult task. Indeed perturbative expansions in SYM assume the 't Hooft coupling 
A to be small, while string perturbation theory makes sense when the string tension 
T = vA/27r is large. A conceptual progress came with the proposal [2] of restricting to 
states (named afterwards BMN states) represented by small closed strings with large 
angular momentum L. Indeed, the energy of these states admits an expansion in the 
small classical parameter X/L 2 , the quantum corrections corresponding to another, 1/L 
expansion. On the SYM side, these string states correspond to composite operators 
containing L local fields. However, this implies that the (for now) accessible semiclassi- 
cal regime of string states is mimicked in the gauge theory by "long" operators, which 
renders a priori the computation of their anomalous dimensions a complex problem. In 
this perspective, it can be understood why a very important progress was realised by 
Minahan and Zarembo while noticing the coincidence between the one-loop expansion 
of the mixing matrix for operators in the M = 4 SYM theory containing L local fields 
and an integrable SO(6), L sites spin chain Hamiltonian For, since then, many 
Bethe Ansatz ideas and techniques were used in order to find anomalous dimensions 
of very long operators of gauge theories, as long as this novel correspondence between 
Integrable Models (IMs) and (S)YM theories was extended to the full PSU(2, 2|4) sym- 
metry jlj, to higher loop expansions 0E1IZ1 and to other gauge theories [KJ EJ Ell HUES- 
For completeness' sake, we should also mention how this IM/SYM relation makes more 
complete and deeper the previous appearance of an integrable system in the barely 
perturbative regime of QCD f3JE]- in particular this link to IMs stimulated an im- 
pressive activity, which allowed many scholars to test the AdS / CFT duality in different 
cases (e.g. for the BMN operators of [2. and for others too). More precisely, the inte- 
grability of the mixing matrix at all orders in perturbation theory was conjectured in 
Ej and then proved for the SU(2) subsector up to three loops in [U]. In this paper, 
the dilatation operator was embedded into a long-range spin chain, the Inozemtsev spin 
chain jT^j. However, the Inozemtsev spin chain at four loops would lose the apparently 
desirable property (in perturbative gauge theory) of BMN scaling and this lack stim- 
ulated the conjecture according to which an alternative long-range spin chain for all 
number of loops may exist [7j. A Bethe Ansatz was also proposed in |7, for all the values 
of the coupling g (cf. below), in order to deal with this otherwise unknown multi-loop 
Hamiltonian. However, the intriguing recent paper [16J has pointed out many reasons 
why this Bethe Ansatz may not furnish the right anomalous dimension at and beyond 
the wrapping order g 2L (note also ^7] for a qualitative interpretation). The same paper 
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has identified in the half filled Hubbard model a short range Hamiltonian conjectured 
to reproduce the mixing matrix. Actually, this identification was explicitly carried out 
up to three loops and it is unclear if it may survive the break-down order g 2L . 

Despite the great amount of work on the subject, the majority of the results found 
up to now (with exception of very few papers like, for instance, [TBI EH EH! ) concerns 
the calculation of only the leading term of the anomalous dimensions of arrays of L 
operators in the limit L — > oo. In the BMN sector, this corresponds to classical ener- 
gies on the string side. Consequently, the correction to this leading term is related to 
quantum fluctuations of the string state energy and is therefore worth studying. In this 
respect and from the Bethe Ansatz point of view the L — > oo limit may be described in 
all its physical quantities in terms of the density of Bethe roots (per quantum numbers), 
provided the latter really tends to a continuous distribution when L = oo (cf. [21] for 
some remarks on this point). In any event, almost as early as the Bethe's invention 
[22] (for the spectrum of the isotropic Heisenberg spin 1/2 chain), a linear integral 
equation constraining the density (for the anti-ferromagnetic ground state) was derived 
and solved [23- Since those early stages the power and versatility of Bethe Ansatz was 
being very much appreciated, at most in condensed matter physics, integrable models 
theory 1 and statistical mechanics (cf. [213 EH] just for some examples). Also, the inte- 
gral equation idea lived a revival since 1964 fZT\ l2~4"| and was extended to excited states 
(cf. the profundus review [28 ) and to the statistical view [21]. Implementing this latter 
in the framework of the relativistic factorised scattering theory, Al. Zamolodchikov for- 
mulated a general and pretty widely applicable idea concerning an exact formula of the 
vacuum scaling function at all size scales, the so-called thermodynamic Bethe Ansatz 

As for finite size effects in quantum integrable systems, the Non-Linear Inte- 
gral Equation (NLIE) description - first introduced in [31] for the conformal (anti- 
ferromagnetic) vacuum and then derived for an off-critical vacuum in [32] by other 
means - turned out to be an efficient tool in order to explore the scaling properties 
of the energy. Since [33], regarding excitations on the vacuum, a number of articles 
was devoted to the analysis of and through a NLIE and mainly follows the route pi- 
oneered by Destri and de Vega [32] (cf- the lectures [31] for an overview). In this 
way (which will be ours too), the NLIE stems directly from the Bethe equations and 
characterises a quantum state by means of a single (or very few) integral equation(s) 
in the complex plane (and possibly some auxiliary algebraic equations). The NLIE 
has been widely studied for integrable models described by trigonometric-type Bethe 
equations: for instance, the 1/2- XXZ spin chain [31], the inhomogeneous 1/2-XXZ and 
sine-Gordon field theory (ground state in [32], excited states in [33]) and the quantum 
(m)KdV/sine-Gordon theory [33] . 

In this paper, we want to propose the Non-Linear Integral Equation idea [3U 1321 

lr The first liaison with the quantum version of integrability, the Yang-Baxter equation, was found 
by the seminal contributions 24 . 
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|33J EU EEj as a tool to compute finite L corrections to the anomalous dimensions 
of (long) operators in M = 4 SYM. In terms of the solution of the NLIE, we can 
indeed write down exact expressions for the observable eigenvalues, as they depend on 
the Bethe roots. Their behaviour for large L may be disentangled analytically and 
numerically. Going into more details, we will concentrate our analysis on the composite 
operator with the highest anomalous dimension. This corresponds in the spin chain 
to the anti-ferromagnetic state, made up of a sea of real Bethe roots. We will also 
study excitations 2 thereof, introduced by the presence of holes. These are the simplest 
possible modifications, as already argued in [HHj, though the anti- ferromagnetic state 
is not the (true) vacuum (with smallest energy or anomalous dimension) of the chain, 
which enjoys a ferromagnetic nature and corresponds, in the gauge theory parlance, to 
the BPS state with all the partons (i.e. the complex scalars) of the same kind. On the 
contrary, it becomes of interest here as its eigenvalue constitutes the upper bound, i.e. 
the largest anomalous dimension: the finiteness of the spectrum is very clear in the spin 
chain and SYM interpretation, although a momentum bound of the string is rather not 
obvious (but semiclassical computation can be trusted in this regime just partially and 
have been started recently jSZ|)- Moreover, it plays the role of the genuine vacuum in 
the large N QCD expansion, £it Insist cit one loop (cf., for instance, JI]). Its interest 
resides also in the fact that the holes excitations will furnish the just smaller anomalous 
dimensions, whose energies are neglected in condensed matter physics since this part 
of the spectrum decouples to infinity from the real spectrum above the ferromagnetic 
vacuum, in the thermodynamic limit. In this respect, we also emphasise that the Neel 
state (| TIT! ••• TO) is n °t an eigenstate in this context. Specifically, we will study 
both the general 50 (6) case (at one loop) and the SU(2) subcase (though this with an 
arbitrary number of loops), providing exact expressions for the anomalous dimensions 
of SYM operators with finite number (L) of fields. In fact, as a starting point for the 
latter we shall use the asymptotic Bethe Ansatz of [Zj, whose reliability is up to order 
g 2L ~ 2 , as already widely stressed. For clarity and simplicity reasons, we will start by 
the exposition of the SU{2) case which, after the proposal by [15], may be interpreted 
as a strong coupling expansion of the Hubbard model at all orders (provided L is large 
enough). We will introduce the techniques in the known example of the Heisenberg 
chain and we will provide original results for the many loop Bethe Ansatz of [7 j , as well 
as for the 50(6) case. 

2 This terminology is borrowed from the cases when the anti- ferromagnetic configuration yields the 
(true) vacuum. 
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2 The SU(2) case: one loop or the Heisenberg chain. 



Let us first consider the SU(2) subsector of the gauge-invariant scalar operators in 
M = 4 SYM field theory. The anomalous dimension of a general composite operator 
containing L scalars is given by 

where A = Ngy M = Sn 2 g 2 is the 't Hooft coupling of the SU (N) super Yang-Mills 
theory and 

M 

E = J2^I (2-2) 

k=l U k + 4 

is the energy of a spin 1/2- XXX chain, i.e. the celebrated Heisenberg spin chain, with 
L sites. Since the pioneering work of Bethe j22] it is well known that the M complex 
numbers (or Bethe roots) must satisfy the equations 

i \ L M 

^4) =l P~ Mfc 7^ J = 1,...,M, (2.3) 

usually named after Bethe as well. In this approach, one set of solutions iden- 
tifies one energy eigenfunction. In the original paper the previous equations are the 
consequence of the imposition of periodicity of the postulated wave function (the fa- 
mous Ansatz), without any clear mention to integrability. This Bethe eigenfunction is 
also (highest weight) eigenstate of the total ^-component spin operator with integer or 
half- integer eigenvalue S = L/2 — M^O. 

Now, we derive a single Non-Linear Integral Equation (NLIE) along the ideas of 
[321 ESj, so that we may have in it a more effective, though equivalent, description 
of Bethe equations. We may need to say that this derivation will have a pedagogical 
purpose in perspective of the multi-loop case of next section, though it will help to 
illustrate the general idea of [321 EH! an d to interpret the results from the gauge theory 
viewpoint. In fact, on the one hand it is just a limiting case of the general Bethe Ansatz 
of the next section, on the other hand similar results are already contained in [3*H I3b|. 

The NLIE will be an equation for the so-called counting function, 

/ i\ M 

Z(x) = L<P (x, - j - fa - u k, 1) > (2-4) 



where the function 



(x, = i In I % + - ) =2 arctan % , £>0, (2.5) 
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is analytic in the strip |Imx| < £ provided the branch of the logarithm is along the 
negative real axis. Then we need a variable that keeps into account the parity of the 
chain in relation with the number of Bethe roots: 

6 = (L-M) mod 2 . (2.6) 

After, by using the simple property 

iln^-iln^=7T, (2.7) 

the Bethe equations can be written in the form 

iL]n^i-J2iln l + Uj - Uk =n(2I 3 + 5-l), j = l,...,M, 

thanks to the introduction of certain integer quantum numbers Ij, or in terms of the 
counting function as 

Z(u j )=ic(2I j + 6-l), j = l,...,M. (2.8) 

The last equations are completely equivalent to the initial Bethe ones (|2.3jh provided 
that Uj enter the counting function by ()2.4j) . 

From now on and only for simplicity reasons we will be considering states char- 
acterised by real roots. This is the formulation proposed in (33] and it can be easily 
extended to deal with arbitrary complex roots. Bearing in mind the limits 

lim 0(ar,O = ±7T, (2.9) 

x— >±oo 

we easily compute the limiting values of the counting function 

lim Z{x) = ±(L - M)?r. 

x— >±oo 

Since Z(x) is an increasing function, the condition ()2.8|) is satisfied by L — M points 
on the real axis, among which there are indeed M Bethe roots. The number of the 
remaining fake 3 solutions (holes) is 

H = L-M-M = L-2M. (2.10) 

Of course, the holes Xh are determined by the same equations as those for real roots, 
but with the complementary set of integer quantum numbers Ih, namely 

Z{x h ) = tt(24 + 5 - 1) , (2.11) 



3 This in the obvious sense that these are not solutions of the initial Bethe equations (|2.3|) . 
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since holes do not satisfy the Bethe equations. Hence, both Bethe roots and holes, 
respectively, enjoy the condition 



exp[iZ(x)] = (—1) , x = Uj,x h . (2.12) 

Now, let f(u) be an analytic function within a strip around the real axis. Thanks to 
(j2.12|) . the sum of its values at all the Bethe roots takes on the expression (33] 

_ M ' r°° Arf 

£/(«*) = -/ ^f'(x-ie)\n[l + (-l) s e iZ ^}- (2.13) 

k=l J 00 

r- 00 ^ JL 
- / ^/'(* + ie) In [1 + (-1)^*^)] - V f(x h ) , 

with e > small enough to keep the integration inside the analyticity strip. If now 
e — ► 0, we may rearrange this expression as 

£/(«*) = " / 5r/'0^)+ (2.i4) 



fe=l 



+ 



.~ 2tt' 

—f\x) Imln [1 + (_l)* e ^(*+")] - ^ /( Xfc ) 

-00 



Upon applying (|2.14jl to the sum over the Bethe roots in the definition of the counting 
function (|2.4|) . we obtain yet a first integral equation for it, 

Z(x) = I4[x,l)-[ ^<P'(x-y,l)Z(y) + 



_ 2tt 

+ — 1mm h + (-l)<V zfe+l0) l + (2.15) 

J-00 
H 

+ y^0(x - 1) . 



h=l 

As usual, we introduce a shorthand 

v(5 JZ(x+iO) 



L(x) = Imln [1 + (_i)V z(x+i0) ] , (2.I6) 
and then take the Fourier transform 4 of all terms in ()2.15|) to obtain 

Z(k) = L0 (k, I) - ±-fr(h, l)Z(k) + -4>'(k, l)L(k) + V e~ ik ^(k, 1) . (2.18) 
v 7 h=l 



l We define the Fourier transform f(k) of a function f(x) as given by 

/>oo 

f(k) = dx e- lkx f{x) . (2.17) 
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This equation can be recast in the more compact form 



H 



Z(k) = F(k) + 2G(k)L(k) + J^ e ~ ikXh H{k) , 

h=l 

where the Fourier transform of the forcing term reads as 

F(k)=L ^ , (2.19) 
1 + i<f>'(k, 1) 

G(k) = ^ilL , (2.20) 
1 + ^0^,1) 

and eventually the holes contribution is (P is the principal value distribution) 

H(k) = _i^lL_ = 2 ^p(\\ G{k) . (2.21) 



that of the kernel as 



^W)=«x^. (2 - 22) 



All these can be easily calculated, once the Fourier transform of the function 

e+x 2 

is explicitly computed as 

0'(A;,O = 27re-«l fc l, (2.23) 

which entails 

/ 1 \ ?7Tf=~^ ~ 1 
F(k) = LP ( — ) , rrrr , Gik) = i"7T . (2.24) 

v ; \k) i(l + e-l fe l) ' V ; l + el fc l V ; 

Upon anti-transforming, we obtain the forcing term 

, * r f 00 dk sinkx r __ L7r _ , 
F(x) = L / r = 2L arctan e — = L gd nx , 



o fc cosh 2 



and besides the kernel 



^ r dk lkx i i d 1 r (i + f ) r (| - 1 ) i d . _ , , 

G(x) = — e lkx ttj = —In —, M- — 7^ N- = — In Six) 

V ' 7_ 00 27T l + el fc l 2nidx r(l-f)T(i + f) 2m dx v ' 

where the expression in terms of Euler's gamma functions, 

r(i + f)r(|-f) 



S{x) 



r(i-f)r(| + f) 



is indeed the scattering factor of the NLIE (cf. EH1 ES| for a justification of this 
name). Finally, from this we can easily gain the hole function in the form 

H(x) = 2ir dyG(y) = -i\nS(x). 



With all these functions at hand, we are in the position to write down the announced 
non-linear integral equation for Z(x), 

Z(x) = F(x)-i^2\xiS(x-x h ) + 2 I dyG(x-y) Imln [l + (-\) s e iZ{y+i0) ] , (2.25) 
h=i 

and we can also check a posteriori that no zero-modes actually entered its derivation. 

The NLIE f!2.25|) together with the holes quantization conditions (j2.11|) is equivalent 
to the Bethe equations ([2.3)1 . 



2.1 The observable eigenvalues. 

Let us now pass on to the computation of the eigenvalues of the observables on states 
containing M real Bethe roots and H holes. We move from 

A r°° dr 

£/(«*) = - / =-f(x)Z(x) + (2.26) 



k=\ 



+ 



—fix) Imln [1 + (_l)* e *^(*-H0)] -J2f(x h ) 

-oo ^ . , 



and then insert into this expression the non-linear integral equation (|2.25j) and re- 
organise the terms as 

M r 00 dr H ( r 00 dr 1 

k=l J 00 71 h=l 71 > 

/°° dr f°° 
—fix) l dy [6(x -y)- G(x - y)] Imln [l + (-l)V z ^ +i0 )] . (2.27) 
■oo ^ J —oo 

This formula gives an exact expression for the eigenvalues of any general observable in 
terms of the solution of the Non-Linear Integral Equation, solution which characterises 
the specific eigenstate. For example, its analogue was exploited in the quantum (m)KdV 
context jHS] to obtain the quintessence of an integrable model, namely the (commuting) 
integrals of motion. Now, we want to use it in order to compute the eigenvalues of the 
energy (anomalous dimension) and of the momentum. 



2.2 The anomalous dimension 

As illustrated in ()2.2|) . to compute the (total) energy, we need to apply formula ()2.27|) 
with the single particle energy 

/(*) = e(x) = -^-j . (2.28) 
x 2 + i 
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Indeed, the first term of the l.h.s. is given by 



00 dr f°° dr f°° dr 1 

^e'(x)F(x) = / ^e(x)F'(x)=L 1 



2n J_ 00 2n J_ OQ 2 x 2 + \ cosh7nr 



00 y 

y 2 + 1 cosh • 



L I dy - ; - _ uny =2Lln2. (2.29) 



The last term reads 



/°° dr r°° 
—e'(x) / dy [5{x -y)- G{x - y)] Imln [l + (-l)^^ )] , 
-00 J -co 

where the ^-convolution is conveniently evaluated in Fourier space (where it becomes 
an ordinary product): 

/°° dn d f°° p ik y 
^ Imln [1 + (-lYe iZ ^ +io n -f / dk— r = 
.00 7T L v ; j dy j_^ 2cosh | 

= f°° dy(4- — !— 1 Imln [l + (-1)V^+^1 . (2.30) 
i-00 \dy coshvn// 

Eventually, we need to compute the two terms of the hole sum (the second term): 

y~N / ^-e'(x) ilnS 1 ^ - - e(a; ft )| . 

For the first of them we may write 

^/t Z" 00 f 00 <7£- 



e\x)ihiS(x-x h )= dxe(x)G(x-x h )= — e lkXh e(k) G(k) . 

2tt y_ oc y_oo 2?r 

This yields, once the second term is expressed by its Fourier transform, 

JL. r°° dk ~ J^L r°° p ikXh 

£ L s ^ 4( *> |GW " 11 = " £ L dk 2^1 ' (2 ' 31) 

where we have used the formula 

e(Jfe) = 2vre^ |fc| , 

particular case of ()2.23|) . and the expression of G(k), ()2.24|) . Eventually, the source 
term may be written as 



y dk Jl y" . (2.32) 

^ ,/_ 00 2 cosh § ^ cosh 7rx h 
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Summing up all the contributions (|2.29l 12.301 12.32)1 , for the eigenvalue of the energy of 
the spin chain, we obtain 

H 

E = 2Lln2-V + 

' COSh TTXh 

h=l rl 

+ r dy( 4- — \ 1 Imln [l + (-l)V^+ i0 )l . (2.33) 

J~oo \dycoshnyJ 

This expression is exact for any L and gives the largest anomalous dimensions of the 
gauge-invariant scalar operators of the SU{2) sub-sector in Af = 4 SYM. The first 
term on the r.h.s. is the known leading term proportional to L; the remaining two 
addends may be expanded in the limit L — > oo to provide respectively 0(1) and 0(1/ L) 
corrections. Analytical expressions up to the order l/L will be given in Section 5. 



2.3 The momentum 

The identification between the anomalous dimension of a gauge-invariant operator and 
the energy of a (spin chain) state needs to be supplemented by the zero momentum 
condition. Therefore, it is necessary to work out, by using the same technology as for 
the energy, the momentum eigenvalue 5 



P 



' M \ 

^p{u k ) J mod 2vr, (2.34) 

\k=l J 



with the single particle momentum, defined as 6 
1 x + - 

p(x) = — In j = 7rsign(x) — 2 arctan 2x = 7rsign(x) — (j) I x, - ) . (2.35) 



x 2 



2 



This relation and the analogous ()2.28|) suggest the interpretation of each root as a par- 
ticle (magnon) exciting the ferromagnetic vacuum and obeying the energy-momentum 
dispersion relation 

eK)=4sin 2 ^, k = l,...,M. 
We remark that p{x) is odd and discontinuous in zero 

p(x) + p(—x) = , lim p(x) = ±n . 

The total momentum may be arranged so as to extract its non-analytic contribution 
P = 7r(M+-M£)-5>U k ,-J =7rM-5>U Jfc) -J , (2.36) 

k=l ^ ' k=l ^ ' 



5 Of course, the identification up to In multiples comes from its definition as a displacement operator 
on a periodic (one-dimensional) lattice. 

6 We need to extend the definition of the sign- function so that sign(0) = 1. 
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with the number of positive or zero real roots and that of negative roots. The 
second equality is obtained, modulo 2tt, by adding 2nM^. Now, we can apply formula 
(I2.27J) to the analytic part of p(x), Pan{x) = ~<\> (x, |). The first term vanishes in that 
the integrand is an odd function 

d t 

=-tf n (x)F(x)=0. (2.37) 
For what concerns the contribution 

— pUx) / dy [6(x -y)- G(x - y)] Imln [l + (-l)V^ )] , 

-oo " J — oo 

we first evaluate the integration on x, 

r°° 

dxp'^(x) [5(x -y)- G(x -y)] = - dk e iky 



oo 



1 + el fc l 



dk 



e iky n 



(2.38) 



2 cosh | cosh ny ' 



in order to obtain 



/°° dr r°° 
— jtjx) / dy [8{x -y)- G(x - y)} Imln [l + (-l)' e **GH-»)] 
■oo ^ J -co 



' — oo 

fOO 



dy — I Im In [l + (-l)V K(2/+i0) l . (2.39) 

cosh ny L J 



Finally, we need to compute the term 

H 



53 {/ w-p'^( x ) i]n3 ( l 

h= l U-oo Z7r 



x-x h )- p m (x h ) \ . (2.40) 



But this expression is the sum of single hole contributions which are minus the primitive 
of ()2.38|) at the value y = Xh- So, each term is given by 

71 . 

arctan smh ny + const. 



cosh ny 



And the integration constant is zero since the single term has to vanish for y = Xh = 
for parity reasons (p an is odd and G even). So (j2.40j) simplifies to 



H { f°° dx 1 
z2\ IT Panfc) HnS{x - x h ) -p an {x h ) > = 22 arctan sinh7rx h . (2.41) 

h=i n > h=i 

All the contributions (|2.371 12.391 12.41|) yield the momentum eigenvalue 

H /"oo -r 

P = 7T M + V (arctan smhnx h ) - / dy — Imln [l + (-l) 5 e iZ(s/+i0) l . (2.42) 

^ J-oo cosh Try 
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The attentive reader will find similar results here and there in [2H EH] • 

As we said at the beginning of this subsection, the condition P = works as a 
constraint for the anomalous dimension (|2.33ll2~Tj) . In particular, the antiferromagnetic 
state simply enjoys 

so it possesses a SYM operator as a counterpart only if L e 4N. 



3 The 577(2) case: multi-loops 

Following the line of the previous Section, we want to establish the NLIE framework 
for the conjectured multi-loop Bethe equations [Zj 



X{ Uj + |) \ L = U j - u k + i 



where we introduced the function 



And as usual the single particle momentum p(uj,X) is such that e ip<yUj ^ L equals the 
l.h.s. of the corresponding Bethe equation (J3.1)) . or explicitly 

1 jC(x+|) 
=-ln— f-. 3.3 

z A (x - |) 

Of course, these equations reproduce the Heisenberg case of Section |21 in the small 
coupling limit A — >• : actually, only the l.h.s. in (j3.1|) . namely the momentum ()3.3|) . has 
changed from the XXX case. Therefore, in the present section we shall systematically 
follow all computations of the previous one. 

In other words, we simply need to change the function (j)(x, |) of ()2.4|) into 



IZ + x) " A 



$ (x,A)=,ln" ; V 4 ^ )2 , (3.4) 

which, contrarily to the momentum ()3.3j) . is continuous in x = 0. Then, the counting 
function may be defined as 

M 

Z(x,A) = L$(x,A)-^0(x-w fc ,l), (3.5) 

k=l 
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so that the Bethe equations read (with certain integer quantum numbers Ij) 

Z(uj, A) = n(2Ij + 8 -1) , (3.6) 
5 = (L - M) mod 2 . 

The asymptotic limits of the counting function are again given by 

lim Z(x, A) = ±(L — M) n 

x— >±oo 

and can be used to fix the number of holes. Indeed, as Z(x, A) is an increasing function, 
the conditions (|3.6|) with generic integers Ij are at most satisfied by L — M points on 
the real axis. This means that the number of holes is 

H = L-2M, (3.7) 

when considering states with real roots only. The position of any hole Xh is fixed by a 
quantisation condition identical to (|3.6jh but with a fake quantum number 

Z(x h , A) = tt(24 + 6 - 1) . (3.8) 

Consequently, both the Bethe roots and the holes, viz. x = Uj, Xh, satisfy the condition 

exp[iZ(x,X)} = (-I)*" 1 . (3.9) 
Again, as e — ► the sum over all the real roots (|2.13|) takes on the form 



k=l 



+ 



E/(«*) = - r ^-f(y)z(y,x)+ (3.io) 

r ^f(y) Imln [l + (-lf e iZ ^} -j^fM ■ 
J -°° 71 h=i 

This expression may be applied to the sum in the counting function ()3.5j) and brings 
Z(x,X) = L${x,\)- [°° ^ ( p'(x-y,l)Z(y,X) + 

+ r —<P'(x-y,l) Imln [l + (-1) ^ e lZ ^ +l0 ^} + (3.11) 

J-oo 7T 
H 

+ ^2<f)(x - X h ,l) . 



h=l 



It is convenient to introduce the usual (cf. (|2.16Jl ) synthetic notation 

L(x,X) = Imln [1 + (_i)V z(x+40 ' A) ] . 
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After x Fourier transforming all the terms and moving the first convolution to the l.h.s., 
we obtain 

H 

Z(k, A) = F(k, A) + 2G{k)L{k, A) + ^ e~ ikXh H{k) . (3.12) 

h=l 

All terms are the same as before in the Heisenberg chain, except the forcing term that 
now depends on A and whose Fourier transform reads 

F{k,X) = L ^ k : X) . (3.13) 

l + i^(M) 

The x Fourier transform of $'(x, A) is given in terms of the Bessel function of the first 
kind J IH 



$'(M) = 2tt e~^J \ ^k 



2tt 



The series expansion Jo(^0kj = (1 — y|^-A + 0(A 2 )) |1U] shows clearly the change with 
respect to the Heisenberg chain. The function 0(x, 1) is unchanged, so we can make 
use of (J2.23|) to arrive at the final expression 

Fl * A) = /. * t 2 col | = ( gd ^ " 16 crf^ + ° (A > ) ' (3 ' 14) 

Inverting the Fourier transforms of ()3.12j) leads to the NLIE valid for this multi-loop 
Bethe equations 7 

H 



Z(x,X) = F(x,\)-iJ2 lnS ( x - x h)+ (3.15) 

h=l 

/oo 
dy G(x - y) Imln [l + (_i)<y zfe+i0 ' A) ] . 
-oo 

Of course, the convolution kernel G and the hole term S are the same as in Section [2J 
what makes the difference is simply the different forcing term F. And besides the 
structure of this NLIE is quite the same as in many other models, except for the specific 
form of the above-computed functions: hence this similarity corroborates straight away 
the effectiveness of our method. 

Therefore, we can still follow the result ()2.27j) on the Heisenberg chain, keeping in 

7 We need to stress anew the absence of an integrable model behind them. 



15 



mind that here the forcing term is given by ([3.14)1 : 

r°° dr 

£/(«*) = -/ ^f'(x)F(x,X)+ (3.16) 
k=i 

—f{x) / cfy - y) - G(x - y)\ Imln [l + (-1) + 

-oo ^ J— oo 

+ ^f'(x) HnS(x-x h )- f(x h ) 

h=i 71 



3.1 Anomalous dimension 

As typical in Bethe Ansatz theory, the energy of the spin chain, and thus the anomalous 
dimension in gauge theory, is given by a sum on all the Bethe roots 

M 

E = Y,e(u 3 ,\), (3.17) 

3=1 

where the (even) single particle energy function equals 

e(x,X) =i\ , X • r ; 1 r- X . (3.18) 

So, we just need to insert this function into ()3.16|) 

r°° dx 

e(w fc , A) = - / —e'(x, X)F{x, A)+ (3.19) 
k=i J -°° 

/OO J /-oo 
— e'(x, A) / cfy [<f(ac - y) - G{x - y)) Im In [l + (_i) V zfe+l0 ' A) ] + 
■oo ^ J —oo 



h=l < J -°<> Z7T 



(x - x h ) - e(x h , A) 

T 

and re-call its Fourier transform 



8vr 2 Ji ( ^k 
e(k, A) =- { T 



y/\ k e ' a' 



In fact, the first contribution reads 



dx e'(x, A) F(x, A) = r p- e(x, A) F'(x, A) 



DO 271 " i-OO 27r 

^ f .^f)^( 2 , 2 -||, + OOT ,. (3,0, 
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Of course, it is the leading term of the anti-ferromagnetic state energy and thus coincides 
with the expression (10) of [TH] (or (17,18) of jH]), where it was interestingly identified 
with the ground state energy of the half-filled Hubbard model. Moreover, the second 
contribution in ()3.19j) may be re-organised by expressing the convolution as an ordinary 
product in the Fourier space: 

/oo f j T poo 
— e'{x, A) / dy [5{x - y) - G(x - y)} Imln [l + (-1) V^ +i0 < A )] = 
-oo ^ J — oo 

oo 

dyet(y, X) Imln [l + (_i)V zfo+i0 ' A) ] , (3.21) 



where we made use of a new function as a Fourier transform of a construct of the Bessel 
function Ji 



d 

dx 



e x {x, A) 
1 



2z Ji ( &k ) e ikx 



dk — _ 

oo v A cosh | 

A d 2 



(3.22) 



cosli7r2 327r 2 dx 2 cosli7ra; 



+ 0(X 2 



The single hole contribution in ()3.19j) is now evaluated by integrating by parts and 
computing the convolution in the Fourier space: 



00 dx 

— e'(x, A) i In S(x — Xh) — e(xh, A) 



2tt 



dx e(x, A) [G(x — Xh) — S(x — Xh)} 



00 dk 

— e lkXh e{k, A) [G(k) - 1] 

-oo 271 



dk2ir e 



ikx h 



Jxl^k 



2tt 



\/A k cosh | 



TT 



A 



cosh 7TXh 32n 2 



d 2 



TT 



+ 0(X 2 ). 



(3.23) 



x=x h 



dx 2 coshvra; 

Eventually, we collect all terms ()3.20l 13.211 13.23|) that form the energy ()3.17|) and obtain 



8tt 



-E 

h=l 



k{e k + 1) 



(3.24) 



dk2ne ikXh _ 
oo v X k cosh TT 



+ 



dy 




2iJ 1 (^k)e ik y^ 

dk x 77 \ 

VA cosh| 



Imln [1 + (-1) V zfe+l0 ' A) ] 



This expression for the anomalous dimension is exact in any regime of L and specif- 
ically the second and third terms provide all the sub- leading corrections when L — > oo, 
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whose expressions are by the same token novel and intriguing. Analytical expressions 
of them up to the order 1/L will be given in Section 5. 

In conclusion, it is worth emphasising that the break-down of the proposal (|3.1|) at 
order g 2L affects our results as a trivial consequence, although it reveals itself unre- 
lated to our method. In other words, the latter should be perfectly applicable to the 
hypothetic correct Bethe Ansatz equations, provided in the typical form, along similar 
steps. 



3.2 Momentum 

We can compute the momentum (|2.34|) by summing the single particle momenta ()3.3|) 
of all the Bethe roots. From ()3.3|) we separate the analytic contribution p an (^fc, A) = 
— $(«fc, A) so that, as in (|2.3fij) . we can write 

M 

P = ttM-^$(u A; ,A). (3.25) 

k=l 

Thus, we only need to apply (|3.16|) to the function p an (x, A). The first contribution 
vanishes 

r°° dx 

-J — y an (x,A)F(x,A) = (3.26) 

as the integrand is the product of an even and an odd function. The second contribution 
to (j3.16|) becomes easily 

/oo J poo 
— p' an (x, A) / dy [5(x - y) - G(x - y)] Imln [l + (-1)V^ ' A )] = 
-oo ^ J — oo 

d JL / dk Imln \1 + (-1)V^+^1 = 

, 7T 2cosh| L v ; j 



oo " ^v"" 2 

oo 



dy F^A) lmln . + ,_ l)5eiZ{yW) , (3 27) 

7T L 

with the appearance of the forcing/momentum term (|3.14j) . The latter also appears in 
the hole contribution to (j3.16J) 

7rP'<m( x i X)ihiS(x -x h ) - p(x h , A) = 

dxp an (x, A) [G(x - x h ) - 8(x - x h )] = ^ ■ (3-28) 

In summary, we have the analogue of (J2.42j) in the present case 

? = ,M + if%A)- f°° dyF^Xl Imln + s e iz( y+i0 ,x)} (3>2g) 
L l i J-oo 71 L - 



h=l 



This completes the general results of the multi-loop scenario and will allow us to extract 
in Section 5 the first finite size corrections analytically (and explicitly). 
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4 The 50(6) scalar sector at one loop: finite size 
results 



We want to illustrate the utility of the NLIE to compute the exact finite size contri- 
butions to the anomalous dimensions and momenta in the SO (6) scalar sector at one 
loop. Therefore, we need to consider a chain of L six-dimensional vectors of the so(6) 
representation. As well known after [12], the Bethe Ansatz diagonalization of all the 
commuting integrals of motion is founded on the following system of coupled Bethe 
equations: 

■ /0 \L Mi . M% . , M s . , 

ttij + 1/2 \ I j" U ltj - U 1)k + l TT U X j - U 2 ,k ~ V2 -|-r Uij - u 3tk -1/2 

u hj -i/2j Ulj _ Ulk _ i 11 Ulj _ U2k + i/2 11 Ulj - u 3<k + i/2 ' 

M 2 ■ Mi . /o 

l = TT U 2 ,j ~ U 2 ,k + I TT" ^2,j ~ U hk - l/2 

Mi 

Ms . • Mi ., 

1 TT M3,j ~ ^3,fc + I TT ~ ^l.fc ~ 

By making use of the function ()2.5|) . we may define three counting functions, i.e. one 
for each group of Bethe equations, 

Mi 



Z x {u) = L(/)(u,l/2)-^2<f>(u-Ui, k ,l) + 

k=l 

M 2 Ms 

+ Y^(P(u-u 2tk , 1/2) + Y^<j)(u-u 3ik ,l/2), 

k=l k=l 

M 2 Mi 

Z 2 (u) = -^0(M-M 2ifc ,l) + ^0( M - Mlifc ,l/2), (4.2) 

k=l k=l 

Ms Mi 

z 3 (u) = - <t>{u - u 3>k , i) + ^(tt - ui,fc, 1/2) , 



^=1 fc=i 

such that the Bethe equations look as 



ZxK,) = ir(2I ltj + 6 X - 1) , 

Z 2 (u2j) = 7r(2/ 2J + 5 2 - 1) , (4.3) 
Z 3 (u 3 j) = vr(2/ 3j + 5 3 - 1) , 
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where Ik,% are integer quantum numbers and the (|2.6jl has been generalised to 

<Ji = (L- M l + M 2 + M 3 ) mod 2 , 

8 2 = (Mi - M 2 ) mod 2 , (4.4) 
§ 3 = (Mi - M 3 ) mod 2 . 

We now consider states of the chain described by real solutions {iik t i} to the Bethe 
equations. For simplicity's sake, we consider the case in which the parities of the integers 
L,Mi,M 2 ,M 3 are such that exp[iZk(uk t i)] = —1 (compare with (|2.12jl ). In addition, 
even if the formalism would allow us to consider states with holes (and generally complex 
roots), here we limit ourselves to states which contain no holes, i.e. such that the points 
u satisfying exp[iZk(u)] = — 1 are exhausted by a (real) solution set {iik,i}: we will be 
extending our results in an incoming publication [13] ■ As before, these requirements 
will constrain the allowed values of Indeed, from the limits ()2.9|) we obtain the 
limiting values 

±tt(L -M 1 + M 2 + M 3 ) , 
±tt(M 1 - M 2 ) , 
±7r(Mi - M 3 ) . 

Imposing the condition that the points u satisfying exp[iZ k (u)] = — 1 are those and 
only those in the solution set {iik,i} furnishes these constraints 

\L - Mi + M 2 + M 3 | = Mi , |Mi-M 2 |=M 2 , |Mi - M 3 | = M 3 , (4.5) 

whose solution, if Mi ^ 0, is 

Mi — L, M 2 = M 3 = ^. (4.6) 

If L e 4N there is one single state with these features and it is the completely anti- 
ferromagnetic state: as discussed in it is the state with maximal energy and zero 
momentum, too. Now, the usual procedure allows us to write a sum over the Bethe 
roots of a function / (analytic within a strip around the real axis) in terms of integrals 
involving the Zs: 

E /(«*,♦) = -/ =-f(x)Zk(x)+ -/'(x)Imln[l + e^^°)] . (4.7) 



lim Z\{x) = 

x— >±oo 

lim Z 2 (x) = 

x— »±oo 

lim Z 3 (x) = 

i— >±oo 
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Upon applying ()4.7j) in the definition of we are on the road to write NLIEs for the 
counting functions: 

Z 1 (x) = Lcf>(x, 1/2) -J ^ 4>'(x - y, 1) Z x {y) + 

+ 2 ^-<p'(x-y, 1) Imln [l + e iZ ^ +i ^} + 
J-oo 2vr 

+ j°° ^ <P'(x - y, 1/2) Z 2 (y) - 2 ^ 0> - y, 1/2) Imln [l + e ^ + ^] + 



JZ 3 (y+iO)] 



+ HP-fix- V, 1/2) Z 3 (y) - 2 r p- 4>\x - y, 1/2) Imln [l + e* 
J-oo 2vr < /_ 00 27r 

Z 2 (x) = -J™ ^^\ x - y ,l)Z 2 (y) + 2j™ ^^\x-yA)Im\n[l + e^y + ^} + 



-oo 



+ f_%^ x ~ v. V2) ^(y) -2 J°° ^-^{x-y, 1/2) Imln [l + e ^ +i0 >] , 

Z 8 (x) = - [°° ^-<f>'(x-y,l)Z 3 (y) + 2 /°° ^ - y, 1) Imln [l + e^^ )] + 
J-oo 2vr J_ 00 2tt 

+ J™ V, 1/2) Zi(y) - 2 ^ <j>'(x - y, 1/2) Imln [l + e **GH-*)] . 

By symmetry considerations (note that M 2 = M 3 = L/2) we can infer that Z 2 (x) = 
Zs(x), so that we have to deal with only two equations. We put again Fourier transforms 
into the game and obtain 

Zi(k) = L$(k,l/2)-^-^'(k,l)Z 1 (k) + 2^-i'(k,l)L 1 (k) + 

271 27T 

+ 2 i- i'(k, 1/2) Z 2 {k) - 4 i- 1/2) L 2 (A;) , 
Z 2 (A;) = - 7 ^0'(M)^(£O + 2 7 ^0'(M)£2(AO + 

Z7T Z7T 



+ i- 1/2) ^i(/c) - 2 i- 1/2) Li(A;) 



and, consequently, 



z l{ k) = ^(M/2) +2 £ 

i + 1 + i) 

+ 2 2 " y V . 7 ; Z 2 (fc) - 4 2?r 7 ; L 2 (fc) , 

i + £tf'(M) i + ^'(M) 

= 2 g^' 1 ) L 2 (k) + ^^ 1/2) Z l( k)- 

i + ^'(M) i + 
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No need to say that the usual short notation 

/oo 
dx e~ ikx Imln [l + e lZ * (x+i0) ] , i = 1, 2 , 
-oo 

has appeared. Thus, clearly the x Fourier transform of the ubiquitous function <f>(x,£), 

#(fc,0 = -27rie-fl*IP^ , (4.8) 



plays a central role to achieve 



Zi(fc) = —L 2ni 6 2 ,,, P ( f ) +2— rrrLi^H 

+ 2 TTFH i2 W- 4 TTFW i2W ' 



Z 2 (k) = 2— rr~L 2 (k)+ C \ u ZAk)-2 e 2 . u LAk). 

zy 1 i_)_ e -|fc| ^ v 1 i_|_ e -|fc| n y i_)_ e -|fe| u 1 

Eventually, these equations can be re-arranged in the clearer manner 

w 1 + e ~\k\ fl\ i _ e -|fc| A 

or, after anti-transforming, in the final NLIEs for the Zs 

/oo 
dyG n (x - y) Imln [l + e iZl{y+i0) ] + 
■oo 

/oo 
dyGi 2 (a:-y)Imlii [\ + e iZ2iy+i0) ] , (4.9) 
■oo 

/oo 
dyG 21 (x - y) Imln [l + e tZl{y+i0) ] + 
-oo 

/oo 
dyG 22 (x-y)lm\n[l + e iZ2{y+i0) ] . (4.10) 
-oo 
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The known functions in the previous equations are explicitly 

j > r f°° „ sin fee cosh | / 

= 2L / dfc — ; = 2L arctan V2 smh 

y re cosh fc v 



, (ire sin rex , , nx 



o 



Gn(x) = G 2 2{x) = - I — coskx — - = (4.11) 

' 2tt cosh re 



o 



i i , i d r(i + f)rQ-f) 

4 cosh 22 2vrz r (l - f ) r (I + f ) 



1,00 dk et 
^12(3^) = 2G 2 i(j:) = — I — - coskx 



7r cosh re 

1 coshf 1 1 : 1 d ln r(l + f)r(|-f) 



V^coshvrx 2cosh7rx 7™ <ix r (| — — ) T (| + — ) 

Summarizing, the equations f)4.9l 14.1(1}) are the Non-Linear Integral Equations describing 
the anti-ferromagnetic state (real solutions without holes to the Bethe equations) of the 
50(6) symmetric chain (vector representation). This state possesses zero momentum 
and the maximal energy. The latter will receive an exact expression - in terms of 
solutions of the NLIEs (|4.9[ 14.10}) - in the next subsection. 

4.1 The anomalous dimension 

The important result of [3] is that the dilatation matrix of scalar operators in M = 4 
SYM at one loop can be mapped to the hamiltonian of an integrable SO (6) symmetric 
chain. In terms of the Bethe roots, its eigenvalue 7 reads as follows: 

A Ml 1 

where E is the chain energy. The maximal eigenvalue (anomalous dimension) is ob- 
tained when considering the solution to the Bethe equations containing real roots and 
no holes. For this configuration, by the same arguments used in the previous sections, 
a sum over the set 1 of Bethe roots can be expressed in terms of integrals involving Z\\ 

J— L r°° /It r°° /It 

£/(«!,<) = -/ pf'(x)Z 1 (x) + 2 p. /'(x)lmln[l + e iz ^)] = 
j =1 J-00 ^ J -00 ^ 



—00 j —00 

dk . . *• f dk 



(27T) 2 ' " ' ^ ' (27T) 



f'(k) Zi(—k) +2 / -—f'^L^-k) 
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Inserting now the NLIE (|4.9|) for Z\ yields 

Upon specializing / to be the single particle energy 

Mi 



E = 2 fl VxT => /(«) = = -^-j , e'(k) = Amk e"^ , (4.13) 
we obtain 



i= i u ii + i m2 + 1 



Z" 00 , 1 + e~ fc z y 00 cosh I . 
E = AL dke~ k - u +2- / rffcfc — — ^ Li(—k) + 

JO 1 + e ^ i-oo cos h « 

+ 2- / dkk—-L 2 (-k). (4.14) 
7r > /_ 00 cosh A; 

The first term in the r.h.s. of (|4.14|) gives the leading contribution when L — ► oo. It 
may be easily evaluated 



4L / dfee- fc 1 + 6 ob =4L / dx^-^-=2L(- + ln2) . (4.15) 
Jo l + e" 2fc 7 1 + x 2 V2 / 1 ; 

We remark that ()4.15|) agrees with (5.10) of [3]. Hence, the maximal energy may be 
re-written as 

. . POO 

E = 2L(- + ln2j+/ dxEi(x)Imln [l + e iZl ( x+uS >] + 

/oo 
^^(^Imln [1 + e iZ2(x+i0) ] , (4.16) 
-oo 

where the last two terms contain the functions 

iM*) = 2V2-^X £ 2 (*) = -7r^L. (4.17) 
ax cosh Tix cosh ^ 

Formula ()4.16|) is an exact expression for the energy in terms of the solution of the 
NLIEs (|4.9| I4.10J1 . When L — > oo, the last two terms provide the 0(1/ L) corrections 
to the anomalous dimension. Analytical computations up to the order 1/L will be the 
topic of the next section. 



5 Analytic calculation of sub-leading order 

It turns out that in all the discussed cases it is even possible to single out the explicit 
sub- leading contribution to the energy as L — > oo. In fact, it is of order 1/L and comes 
out in a rather standard way by following the strategy of the 'derivative lemma' f32| . 
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Also, we need to mention that higher order corrections might be extracted explicitly 
in the framework of NLIEs; for Klumper was able to compute the first logarithmic 
corrections (to the 1/L term) in the spin 1/2-XXX Heisenberg chain starting by a 
NLIE similar to ours, although using some numerical insights [36.. (see also jUj and 
references therein for a comparison with computations in a field theoretic framework). 
For the time being, we are not interested in discussing the analytic derivation of these 
contributions, but will return to them in a subsequent paper |43j . Nonetheless, their 
presence too is expected in the multi-loop SU(2) and £0(6) cases, as motivated in 
Section 6, where we will present our numerical results. 



XXX Heisenberg chain. Let us first consider the Heisenberg chain and the 
excitations over the anti-ferromagnetic state described by holes. The contributions we 
want to evaluate are contained in the integration term of ()2.33j) : 

sinh Try 



AE(L) 



dy (-tt) 



cosh Try 



Imln [1 + (-l)V zfe+l0) ] 



d y (_„) J^llIlL Imln [ X + ( _ 1)V Z (j/ +iO)] + 



oo 
oo 



+ / dy(-n[ 



cosh 2 Try 
sinh Tiy 



(5.1) 



Imln[l + (-l)V zfe+i0) ] 



cosh Try 

In order to single out the order 1/L contributions, we perform different changes of 
variables in each of the integrals in (j5.1j) , 

n ln2 L r „ ln2L 

y = , for y < , y — 9 H , for y > , 



7T 

and then we let L — > oo. In this limit we have 

'•OC 



7T 



AE(L) 



TT 

L 

+ o(l/L) 



d9U ne Imln 



1 + e 



iZ~(0+iO) 



-e~ nB Imln 



1 + e 



iZ+(8+iO) 



+ 

(5.2) 



where the usual symbol o(l/L) is used to indicate terms that vanish faster than 1/L 

lim o(l/L)L = 0. (5.3) 

L^oo 

The new functions 



Z T (6) = lim 



Z 6 T 



In 2L 



7T 



7T 



±-(H + L) 



satisfy the kink NLIEs: 

Z~{6) = e n9 + 2 



d& G{0-&) Imln 



1 + e 



iZ-(6'+iO) 



(5.4) 



Z + (9) = -e~ n9 + 2 dff G(e-e') Imln 



1 + e 



iz+(e'+io) 
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We notice that the holes contribution in the NLIE ()2.25|) has become a constant and 
has been reabsorbed in a redefinition of the counting function. 

In summary, the integral term in ()5.2|) gives the contribution of order 1/L to the 
energy (as L — > oo). This term can be exactly computed by using the derivative lemma 
based on the dilogarithmic function (for an enunciation of this lemma see for instance 
Section 7 of [H2])- The consequent result for AE(L) is 

AE(L) = ^ + o(l/L) (5.5) 

and it does not depend on the holes, as long as their position remains finite for large L. 

We can now evaluate the holes contribution to the energy ()2.33j) . From 1)2.11)1 we 
deduce the leading behaviour 

x h ~ 7 , (5.6) 

namely holes accumulate towards the point x — 0. This leads to the following contri- 
bution, 

- t ^ = -» * + W X> + S - If + 0(1/L-) , (5.7) 

h=l h=l 

from which we conclude that hole excitations over the anti-ferromagnetic state do not 
contribute to 1/L terms. In summary, the energy in given by 

7T 2 

E = 2L In 2 - ttH + — + oil L) , (5.8) 
6L 

where the symbol o(l/L) is defined in ()5.3j) . 

The SU(2) multi-loop chain. In this case we will evaluate explicitly the energy 
(|3.24|) up to the order 1/L (in the limit L — > oo). Let us start by the term, 

/oo 
dy ei (y, A) Imln [l + (_i) V z(rH0 ' A) ] , (5.9) 
-oo 

where ei(x, A) has been defined in ()3.22|) . As before, we split the integral ()5.9j) into two 

parts, 

r° 

AE(L, A) = dy ei (y, A) Imln [l + (_i)V z(y+i0 ' A) ] + 

J — oo 
/■oo 

+ I dy e x {y, A) Imln [l + (_i)V^ +i0 ' A )] , 
Jo 



we perform the change of variables, 



y = e-^^, fory<0, y = e+ 1 -^^, foryXD, (5.10) 

7T 7T 
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and then we let L — > oo. In this limit the two integrals can be computed by using the 
residue method. The first (second) one is evaluated after choosing a contour closing in 
the lower (upper) complex y-plane. The poles of the integrand lie on the imaginary axis, 
k = ±iii(2n + 1) , n > and give contributions proportional to L~ 2n ~ l . Restricting to 
the leading 1/L contribution in the limit L — > oo, we can write 



dd e ny 1mm 



1 + e iZ-(9+iO,\) 



(5.11) 



A-Ki ( ^fk 
j=Ji —i — 



dde'^ Imln 



1 + e iz+(e+io,\) 



+ o(l/L) 



where o(l/L) satisfies (J5.3)) and, as before, we have defined the new functions 

ln2L 



Z^{6,\) = lim 

L^oo 



Z 9t 
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-,\\±-(H + L) 

71 I 



(5.12) 



The equations satisfied by Z T (6, A) are obtained starting from (J3.15|) . performing the 
shifts appearing in their definition ()5.12j) . then evaluating the leading contribution of 
holes and the forcing term when L — > oo. The holes term gives, as before, ^ttH/2. On 
the other hand, the forcing term can be evaluated by using the residue technique in a 
fashion similar to the previous energy kernel calculation. Its contribution is 



± Jo (i^ ] + 0(1/L) 



(5.13) 



the first term coming from the residue in k = 0, the second from the residues in k = ±i7T 
of the integrand of ()3.14|) . It follows that the equations satisfied by Z^(9, A) take the 
form 



Z-(6,\) = J ( ^ + 2 / d9' G{9-6') Imln 



1 + e 



iZ-{6'+iO,X) 



(5.14) 



Z+(9,\) 




60 G(6-6') Imln 



1 + e 



iZ+(0'+iO,\) 



We remark that we are still able to apply the derivative lemma in order to compute the 
1/L contributions of (|5.11|) . the only difference with the Heisenberg chain being that 
now the functions involved are dressed with Bessel functions. The result is 



AE(L, A) 



At J i(-^)tt 2 
L ^ Jo (if) 6 



o(l/L). 



(5.15) 
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We now consider the holes contribution in (|3.24|) : 



H ,00 J, (&k) 

Y dk2ire ikXh ^ '—. (5.16) 

^ 1 v/A^coshl 



h=l J -°° v " 2 

From (J3.8)) we deduce again the leading behaviour 

rr h ~ 7 ■ (5.17) 

It follows that the holes contribution to the energy is 

where £^ 2 z are the Euler numbers. 

Therefore, summing up ()5. 151 15.18J1 with the leading contribution to (|3.24|) propor- 
tional to L, we get that the energy of the multi-loop chain in the limit L — > 00 behaves 
as follows 

r^O^H^) (-a)' 



E = —!=L\ dk — ^—4 \ - Htt > tt~~t: ; - LEa + 

v/A 7o A;(e fc + 1) ^ 2 *Z!(Z + 1)! 1 2/1 



4i J i(-^)vr 2 . 1/r . 

+ ^Tpf^ +0(1/L) - 

The 50(6) symmetric chain. In an analogous way, one can estimate exactly the 
coefficient of the 1/L correction to the energy (|4.16|) of the 50(6) chain (as L — > 00). 
This correction is contained in the two integrals of ()4.16|) 

AE(L) = J2 d V E i(y) Imln I 1 + e lZ > {y+i0) ] = 
i=i 
2 ,0 

= X) / ^-^(^ Imln [ X + e lZ ^ y+l0) ] + (5.19) 

2 ^oo 

+ V / dyEi(y) Imln [l + e 4Z ^ +i0 )] . 
i=i Jo 

As in the 577(2) case, we perform different changes of variables according to the region 
of integration, 

2 2 
y = In 2L , for « < , y = + - In 2L , for y > , 

7T 7T 
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and then we let L — > oo. In this limit we get 

AE(L) = 

+ e 2 1mm 

2 l J 

rr / , | 1 



VT 1 jr9 

— / afc' < — = e 2 Im In 



(5.20) 



+- 



dO 



V2 



e 2 Im In 



1 _|_ e ^r(e+iO) 
1 + e iz+(e+io) 



1 + e 



iZ~(0+iO) 



+ 



— e 2 Im In 



1 + e 



+ 



+o(l/L) , 

where o(l/L) indicates terms that vanish as (|5.3|) . The new functions 



Zj(0) = lim 

L— >oo 



, Z*(0)=lim 

L— >oo 



Zi ^0 + — ln2Lj ±Ltt 
satisfy the /cznA; NLIEs: 

/oo 
d0' GniO-O 1 ) Imln 
■oo 

/■oo 

+ 2 / d0' G 12 (0-0') Imln 



Z 2 =F - In 2L ± L 



7T 



7T 



1 + e 



\zf(e'+io) 



1 + jZ^QV+IO) 



(5.21; 



ZJ(0) = ±e^+2 dO' G 21 (0-0') Imln 



1 + e 



i£f (0+iO) 



+ 



+ 2 / d0' G 22 (0-0') Imln 



1 + e ^ 2 T (e'+io) 



Now, it happens that the two contributions of order 1/L in (|5.20J) can be exactly 
computed by generalizing the derivative lemma to the SO (6) with two 

coupled NLIEs. These two contributions are equal and their sum gives the order 1/L 
contribution to the energy: 

(5.22) 



7T 



AE(L) = — + o(l/ L) . 

We conclude that in the limit L — > oo the energy (|4.1fi|) of the anti-ferromagnetic state 
of the 5*0(6) chain is given by 

E = 2L{^- + \n2)+^ L + o{l/L). (5.23) 

This finite size correction induces to think of a c = 3 two-dimensional conformal field 
theory. 



6 Numerical analysis 

In the previous section we have computed the explicit contribution to the energy up to 
order 1/L for growing L. However, the NLIE formulation of the Bethe Ansatz equations 
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allows us to work out interesting and precise numerical computations as well (cf. [BIH ISH] 
as first numerical results within this approach). The latter can be used, for instance, to 
confirm and to improve analytical results. In this spirit, we have performed numerical 
calculations in order to study contributions to the energy of orders equal to and smaller 
than 1/L, when L is very large. 

We show here few examples of our numerical solutions of the equations. The most 
natural method is to solve them by iterations [3B]- Even if the obtained equations are 
correct for all numbers of holes and lengths of the chain, they are particularly effective 
if one considers states with a small number of holes in a 'Fermi-Dirac sea' of real roots. 

Heisenberg chain. We consider a zero momentum state with four holes quantised 

-0.2 | 1 1 1 1 1 



-0.25 - 




Figure 1: Plot of AE(L)L — ^- versus L for the state I = (—1, 0, 1, 2) of the Heisenberg 
spin chain. 

by I — (—1, 0, 1, 2) (see also Fig. |21 for a prototypical example of the behaviour of 
the counting function) and we follow the evolution of E with L. The goal is to show 
the order of contribution of the various terms in (J2.33)) when L is large. The leading 
contribution is explicit: 2Lln2. From the discussion in the previous section, we know 
that the holes depending term gives a contribution — Htc + 0(1/L 2 ). The contribution 
in 1/L is contained in the integral term (j5.1|) that behaves as 

AE(L) = ^ + o(l/ L) , (6.1) 

where o(l/L) indicates corrections that vanish faster than 1/L. While trying to obtain 
some insight on them we extended our analysis to chains of up to five millions sites 
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2 

because we observed that the quantity AE(L) L — ^- is extremely slow to vanish at 
growing L, as clarified by Fig. ^ 

This behaviour would be perfectly consistent with the presence of logarithmic terms 
(inside o(l/L)). Further, we expect the first one of them to take the following form 



7T 



AE(L)L- — 

6 



Cl 



ln^L 



+ . . . 



(6.2) 



and we have found a good agreement with the data coming from the numerical solution 
of the NLIE. However, we refrain from giving any estimate of the constant ci, because 
numerical computations for very large chains are technically difficult and are affected 
by growing numerical errors. We expect to produce more precise data in the future to 
better understand the additional terms in (|6.1|) . 




Figure 2: Plot of the counting function Z(x) versus x for the Heisenberg spin chain 
with L — 12 sites. The position of the four holes quantised by I = (—1, 0, 1, 2) is 
indicated by the small crosses. 



The finite size corrections to the anti-ferromagnetic vacuum of the Heisenberg chain 
have been extensively studied in condensed matter literature EQ. Taking into 
account the first logarithmic corrections we have 



E = 2L ln2 



7T 

QL 



1 + 



3 1 



81n J L 



lnlnL 

+ k 2 ^^- 4 



ln 4 L ln 4 L 



+ . . . 



(6.3) 



where the numerical values of the constants k 2 and k% have been obtained by a fit of 
the numerical data and have been found to be in agreement with those of Karbach 
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and Mutter [13]. We remark that logarithmic corrections to the energy of the anti- 
ferromagnetic state are smaller than those for states containing holes. This seems to 
be a general property of this state. 



SU(2) at many-loops. We can perform similar investigations on AE(L, A) in (|5.9|) 
when L is large. Even in this case, (with or without holes) we immediately see that 
o(l/ L)L vanishes very slowly with L and this suggests the presence of logarithmic 
corrections. As pointed out in [HB] this logarithmic behaviour is related to the decay 
at infinity of the kernel as a power. Since the multi-loop Bethe equations of [Zj provide 
the same kernel as the Heisenberg chain, the presence of such logarithmic corrections 
is somehow expected. 

We used numerical data, obtained for the same four hole state we used in the 
Heisenberg case and for A = 50, to test the agreement with the following guess for the 
logarithmic corrections 



AE(L, A) L - 0.7843670037. 



bi 



\n 2 L 



+ ... . 



(6.4) 



(on the left hand side we have provided the numerical value of the coefficient of 1/L in 
(|5.15|) for A = 50). As in the Heisenberg case, we refrain from giving any estimate of 
the constant b±, because of the growing numerical errors that affect computations for 
very large chains and we postpone them to a forthcoming publication 
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Figure 3: Plot of AE(L, A) L versus L for the state I = (—1 , , 1 , 2) of the many-loops 
spin chain with A = 50. 



However, in close analogy with the Heisenberg chain, the proposed functional form 
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(j6.4j) for the logarithmic corrections was found to be in good agreement with the nu- 
merical solution of the NLIE. In Fig. |3] it is possible to find the plot of the data for 
AE(L, A) L which have been used for our analysis. 

Finally, we guessed that the finite size corrections for E(L, A) in the case of the 
anti-ferromagnetic state will have the same structure of those corresponding to the 
Heisenberg chain 



VX Jo Ke k + 1) LV\ J ' 



o(l/L) = T^h + ... • (6.5) 



with 

c(A) 
L In 3 L 

Again, in the case A = 50, we obtained a good agreement between the previous expres- 
sion and the numerical data. 



The SO (6) chain. The same procedure used for the SU{2) case applies to the 
SO (6) chain, the main difference being that now we need to solve two coupled equations. 
This does not imply any additional difficulty in their solution, but simply increases the 
computational time, and hence we prefer to limit our analysis to the anti-ferromagnetic 
state for now. As before, we investigate AE(L) L for large L and compare numerical 
results (coming from (|4.16Jl ) with the analytical evaluation (|5.22|) . according to which 

AE(L)L- — = Lo(l/L). (6.6) 

Thanks to the various insights from the above-discussed cases, we may for now limit 
our analysis to chains with up to L = 1200 sites: we intend to improve the details 
in an incoming paper along with analytic supports. Nevertheless, although not 
completely definitive, the leading behaviour of Lo(l/L) seems to repeat that of the 
Heisenberg ground state, namely l/ln 3 L. The corresponding picture is in Fig. HJ 



7 Conclusive remarks 

Briefly said, the novelty of this work may be summarised in the derivation of the NLIEs 
and their application to computing anomalous dimensions and momenta in M = 4 
SYM theory when the number of operators, L, is finite. More specifically, our attention 
has been focussed on the SU(2) scalar sub-sector for many loops and the more general 
5*0(6) scalar sector for one loop. And for what concerns the SU(2) case we started by 
the conjectured Bethe equations of [Zj, which are believed to be trustable at least up 
to order g 2L ~ 2 . Actually, we have also re-derived, following our own route, the know 
formulae for the isotropic spin 1/2 XXX Heisenberg chain as finite L results for the one 
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Figure 4: Plot of AE(L) L — ^ versus L for the anti-ferromagnetic state of the SO(6) 
spin chain. 

loop correction in the SU(2) scalar sub-sector: this effort was at least useful to warm 
up and check the entire machinery. 

These NLIEs, equipped with the quantisation conditions for the holes, are totally 
equivalent to the Bethe equations we started with. But of course they are much more 
effective for both the analytic and numeric computations. For simplicity's sake, we 
limited our analysis to the case of real roots, although the introduced formalism can 
be easily extended to allow for the possible presence of complex roots, as initiated by 
[SHI EHI - Of course, because of its intrinsic nature this NLIE idea is of easier applicability 
in presence of a very large number of real roots, a limited number of holes and possibly 
of complex roots. In terms of the spin chain this corresponds to focus the attention on 
the anti-ferromagnetic state and excitations over it. However, the anti-ferromagnetic 
configuration may receive particular interest in SYM theory as that with the largest 
anomalous dimension; so, in this respect it is antipodal to the ferromagnetic vacuum. 

Nonetheless, the present formulation gives rise to compact and exact expressions 
for the 'observables' valid for any L. For instance, the anomalous dimension assumes 
a form which, in principle, could be exactly computed after solving the NLIE for the 
counting function Z(x) and after fixing the hole positions, Xh, by the quantization 
conditions. In this respect, the analytical results of Section 5 have given explicitly the 
anomalous dimension up to the order 1/L (in the limit L — > oo). In addition, the 
numerical work of Section 6 has also shown the next-to-leading logarithmic dependence 
(on L). In principle, the same ingredients might also be exploited to obtain all the other 
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conserved charges underlying integr ability [33]. Therefore, a comparison with the string 
theory integrability, first disclosed by [46J, would be highly instructive and desirable. 

On the other hand, we would like to remark that the method used here is quite 
flexible and can be applied to various integrable models. Among them it is important 
to mention the Hubbard model for its recent relevance in the computation of anomalous 
dimensions of SYM m this respect, it is not clear whether the Hubbard chain might 
be the ultimate model in the SYM/IM relation since there seems to be no trace of such 
a model (or its symmetries) on the string side (not to say on the SYM side). 

Finally, it will be interesting to apply the present approach to the case of large- ./V" 
QCD [TT], where the integrable anti-ferromagnetic spin-1 XXX chain appears in the 
computation of anomalous dimensions. 
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